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Abstract 


For the 0(A) field theory with A<3? 4 self-coupling, we construct the two- 
particle-irreducible (2PI), closed-time-path (CTP) effective action in a general 
curved spacetime. From this we derive a set of coupled equations for the mean 
field and the variance. They are useful for studying the nonperturbative, 
nonequilibrium dynamics of a quantum field when full back reactions of the 
quantum field on the curved spacetime, as well as the fluctuations on the mean 
field, are required. Applications to phase transitions in the early Universe such 
as at the Planck scale or in the reheating phase of chaotic inflation are under 
investigation. 
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I. INTRODUCTION 


One major direction of research on quantum field theory in curved spacetime ]T]-Q since 
the 1980s has been the application of interacting quantum fields to the consideration of 
symmetry breaking and phase transitions in the early Universe, from the Planck to the 
grand unified energy scales [|4HT2|. In a series of work, Hu, O’Connor, Shen, Sinha, and 
Stylianopoulos [fHil [JT]] systematically investigated the effect of spacetime curvature, dynam¬ 
ics, and finite temperature in causing a symmetry restoration of interacting quantum fields 
in curved spacetime. In general one wants to see how quantum fluctuations 99 around a 
mean field (j) change as a function of these parameters. For this purpose, the two-particlc- 
irreducible (2PI) effective action was constructed for an iV-component scalar O(iV) model 
with quartic interaction |T~^.[T0|.p~D| . Hu and O’Connor [TIJ found that the spectrum of the 


small-fluctuation operator contains interesting information concerning how infrared behavior 
of the system depends on the geometry and topology. The equation for </> containing con¬ 
tributions from the variance of the fluctuation held ( 99 s ) depicts how the mean held evolves 
in time. This program explored two of the three essential elements of an investigation of a 
phase transition [|16j, the geometry and topology and the held theory and infrared behavior 
aspects, but not the nonequilibrium statistical-mechanical aspect. 

For this and other reasons, Calzetta and Hu [[?2| started exploring the closed-time-path 
(CTP) or Schwinger-Keldysh formalism |23|-f26|j, which is formulated with an “in-in” bound¬ 
ary condition. Because the CTP effective action produces a real and causal equation of 
motion [^7|,|8|, it is well suited for particle production back-reaction problems j[2fj| pTT| . Use 
of the CTP formalism in conjunction with the 2PI effective action and the Wigner 


function enabled Calzetta and Hu to construct a quantum kinetic held theory (in hat 
spacetime), deriving the Boltzmann held equation from hrst principles [[34]. The necessary 
ingredients were then in place for an analysis of nonequilibrium phase transitions |35|. In 
recent years these tools (CTP, 2PI) have indeed been applied to the problems of heavy-ion 
collisions, pair production in strong electric helds |36|], disoriented chiral condensates |37| , |38| , 
and reheating in inflationary cosmology |59|. However, none of these recent works has in¬ 


cluded curved spacetime effects in a self-consistent manner, where the spacetime governs the 
evolution of a quantum held and is, in turn, governed by the quantum held dynamics. This 
is especially important for Planck scale processes involving quantum fluctuations with back 
reaction, such as particle creation |4(J, galaxy formation |4Tp, preheating, and thermalization 


in chaotic inflation 42,43 


With this paper we return to the problems begun by Calzetta, Hu and O’Connor a decade 
ago. We wish to derive the coupled equations for the evolution of the mean held and the 
variance for the O(iV) model in curved spacetime, which should provide a solid and versatile 
platform for studies of phase transitions in the early Universe. The hrst order of business 
is to construct the CTP-2PI effective action in a general curved spacetime. The evolution 
equations are derived from it. We must also deal with the divergences arising in it. From 
the vantage point of the correlation hierarchy (and the associated master effective action) as 
applied to a nonequilibrium quantum held |44|, there is a priori no reason why one should 
stop at the 2PI effective action. Indeed, the 2PI effective action corresponds to a further 
approximation from the two-loop truncation of the master effective action constructed from 
the full Schwinger-Dyson hierarchy [|45|,[b|] . For problems where the mean held and the two- 
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point function give an adequate description (which is not the case near the critical point, 
where one has to be careful), the CTP-2PI effective action is sufficient. In particular, the 
2PI effective action contains the commonly used large- N, time-dependent Hartree-Fock, and 
one-loop approximations. 

The O (N) model has been usefully applied to a great variety of problems in field theory 

46fl . The O (N) field theory has the advantage that it affords 
[fTipbl], which yields nonperturbative evolution equations in the 


and statistical mechanics 
use of the 1/N expansion 


regime of strong mean field (it yields local, coupled dynamical equations for the mean 
held and the mode functions of the fluctuation held). Recently it has been applied to 
problems of nonequilibrium phase transitions |47|—f49| . In the preheating problem studied in 
the following paper ]50|], we shall see that this is particularly important for chaotic inflation 
scenarios |d]], in which the inhaton mean-held amplitude can be on the order of the Planck 
mass at the end of the slow roll period [|52l, |53|. The 1/N expansion has many attractive 
features, as it is known to preserve the Ward identities of the O (IV) model |53|] and to yield 
a covariantly conserved energy-momentum tensor |]5fJ. Furthermore, in the limit of large 


N, the quantum effective action for the matter helds can be interpreted as a leading-order 
term in the expansion of the full (matter plus gravity) quantum effective action |55| 


Mazzitclli and Paz J72J have studied the A<h 4 and O(IV) held theories in a general curved 
spacetime in the Gaussian and large -N approximations, respectively. Their approach differs 
from ours in that it is based on a Gaussian factorization which does not permit system¬ 
atic improvement either in the loop expansion or in the 1/IV approximation. In contrast, 
our approach is based on a elosed-time-path formulation of the correlation dynamics, and 
the evolution equation we obtain for the two-point function contains a two-loop radiative 
damping contribution which is not present in the large-IV approximation. At leading order 
in the large-IV approximation, our results agree with theirs, so that their renormalization 
counterterms can be directly applied to the mean held and gap equations derived here. 

This paper is organized as follows. Secs. |Tfj and H present self-contained summaries 
of the two essential theoretical methodologies employed in this study, the closed-time-path 
formalism and the two-particlc-irreducible effective action. The adaptation of these tools to 
the quantum dynamics of a A<F 4 held theory in curved spacetime is presented in Sec. m 
The O(IV) scalar held theory is treated in Sec. 0. 

Throughout this paper we use units in which c = 1. Planck’s constant h is shown 
explicitly (i.e., not set equal to 1) except in those sections where noted. In these units, 
Newton’s constant is G = hM P ~ 2 , where M P is the Planck mass. We work with a four- 
dimensional spacetime manifold, and follow the sign conventions^ of Birrell and Davies |2| 
for the metric tensor g IUJ , the Riemann curvature tensor Rp Ua p, and the Einstein tensor Gp V . 
We use greek letters to denote spacetime indices. The beginning latin letters a, 6, c, d, e, / 
are used as time branch indices (see Sec. 0 ), and the middle latin letters i,j,k,l,m,n are 
used as indices in the O(IV) space (see Sec. 0)- The Einstein summation convention over 
repeated indices is employed. Covariant differentiation is denoted with a nabla or a 


1 In the classification scheme of Misner, Thorne and Wheeler [Q, the sign convention of Birrell 
and Davies Q is classified as (+,+,+). 


3 





























semicolon. 


II. SCHWINGER-KELDYSH FORMALISM 


The Schwinger-Keldysh or “closed-time-path” (CTP) formalism is a powerful method 
for deriving real and causal evolution equations for expectation values of quantum operators 
for nonequilibrium fields, i.e., for quantum systems where the density matrix p and the 
Hamiltonian H do not commute, [H, p] 7 ^ 0. This can occur, for example, in a field theory 
quantized on a dynamical background spacetime and also in an interacting held theory with 
nonequilibrium initial conditions. The methods discussed here are well suited to studying 
the dynamics of an open quantum system. Excellent reviews of the Schwinger-Keldysh 
method are Zhou et al. ]26|] as applied to nonequilibrium quantum held theory and Calzetta 
and Hu |22j as applied to back reaction in semiclassical gravity. In this section we briehy 
review the Schwinger-Keldysh method in the context of an interacting scalar held theory in 
Minkowski space, with vacuum boundary conditions. 

Consider a scalar held $ in Minkowski space with a A $ 4 self-interaction. Studying the 
semiclassical properties of the theory consists of taking the degrees of freedom to be the 
classical held cj) and huctuations ip about the classical held conhguration. The equation of 
motion for small oscillations of (j) about the stable quantum-corrected equilibrium conhg¬ 
uration is obtained via a variational principle from the effective action T[0] fl58|. In the 
conventional Schwinger-DeWitt or “in-out” approach [lj.|57|], one couples a c-number source 
J (which is a function on M 4 ) to the held </> and computes the vacuum persistence amplitude 
in the presence of the source J. This amplitude has a path integral representation 


Z[J\ = exp CjW[J] 


= / D(j)ex p 


S F [(j)\ + / d 4: xJ(x)(j)(x] 


( 2 , 1 ) 


where the functional integral is a sum over classical histories of the (j) held which are pure 
negative frequency [i.e., all spatial Fourier modes of </> have a time dependence like exp (iut), 
to > 0 ] in the asymptotic past and pure positive frequency [~ exp(— iut)\ in the asymptotic 
future . 0 

I 11 a general nonequilibrium setting, such as in a curved or dynamical spacetime or 
when (j) is time dependent, the notion of positive frequency in the asymptotic past is in 
general different from that in the asymptotic future. Hence, the “in” vacuum state 10, in) 
dehned at x° = —00 and the “out” vacuum state | 0 ,out) defined at x° = 00 are not 
necessarily equivalent. The generating functional Z[J] dehned in Eq. (2.1) is then the 
vacuum persistence amplitude 


( 0 , out| 0 , in)j = ( 0 ,out|Texp ( — / d A xJ(x)$ H (x) ] |0, in), 


( 2 . 2 ) 


2 It is noted that these boundary conditions on the functional integral are equivalent (up to an 
overall normalization) to adding a small imaginary term —iecj) 2 to the classical action, where e > 0. 
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where $ h ((e) is the Heisenberg held operator for the theory without the source J. This 
amplitude is in general complex. It follows that the classical held obtained by functional 
differentiation of — ihlnZ[J] is the matrix element (0, out|$ H |0, in) which will in general be 
complex. In addition, the dependence of 0j = SW/6J on J will not, in general, be causal 
p8| , [27[| . In curved spacetime, the energy-momentum tensor (T0,) is obtained by functional 
differentiation of W with respect to , which at one loop yields a complex matrix element 
of X0,(<1> H ) between the “in” and “out” vacua, where <h H is the Heisenberg held operator and 
X0,(0) is the classical energy-momentum tensor for the held [|l|. 

In the closed-time-path formalism, real and causal dynamics for 0 can be obtained, as 
well as the expectation value of the energy-momentum tensor. Let x° = x ° be far to the 
future of any dynamics we wish to study. It is not necessary to assume that A = 0 or that 
the Hamiltonian is time independent at x° = x®. As in the previous “in-out” approach, 
suppose we wish to compute the quantum-corrected equation governing the classical held 0. 
Let M = {(:r 0 , T)| — oo < x° < be the portion of Minkowski space to the past of time 
x*. We start by defining a new manifold as a quotient space 

M = (Mx {+,-})/~, (2.3) 

where ~ is an equivalence relation defined by the following rules 

(x, +) ~ (x', +) iff x = x 1 

(x, —) ~ ( x', —) iff x = x' (2.4) 

(x, +) ~ (x', —) iff x = x' and x° = x®. 

The manifold JA is orientable, provided we reverse the sign of the volume form between the 
+ and — pieces of the manifold. It is then straightforward to generalize the usual effective 
action construction to the new manifold JA. With the volume form on JA, we can generalize 
the classical action S F to JA , 

S F [0+, 0-] = S F [<P+] ~ S F [0-], (2.5) 

where S F [(p] is the classical action on M, and <p + and 0_ denote the 0 held on the + and — 
branches of JA, respectively. The spacetime integrations in the right-hand side of Eq. (|2.5|) 
are understood to be over M. In order for <fi± to be a function on JA, we must have 


0+0*01*0 = 0-0*013.0. (2.6) 

The generating functional of vacuum n-point functions (i.e., expectation values in the 10, in) 
vacuum) for this theory is then dehned by 


Z[J + , J_] = j D(j) + D(j)_exp 
J ctp 


- ( 5 F [0+,0_] + / d x(J+0 + - J_ 


'M 


(2.7) 


where J + and J_ are c-number sources on the + and — branches of JA, respectively. The 
designation “ctp” indicates that the functional integrals in Eq. ( |2.7| ) are over all held con¬ 
figurations (0_|_, 0—) such that (i) 0 + = 0_ at the x° = x ° hypersurface and (ii) 0 + (0_) 
consists of pure negative (positive) frequency modes at x° = — oo. It is not necessary for the 
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normal derivatives of 0 + and 0_ to be equal at x° = x*. Because the theory is free in the 
asymptotic past, a positive frequency mode[] is a solution to the spatial-Fourier transformed 
Euler-Lagrange equation for 0 whose asymptotic behavior at x° = — oo is exp(—Lux 0 ), for 
u > 0 . 

The generating functional for connected diagrams is then defined by 

W[J+, J_] = -ihlnZ[J + , J_], (2.8) 


Classical helds on both + and — branches are then defined as 



j, 5W[J+,J-] 

5J b (x) 


(2.9) 


where a, b are time branch indices with index set {+,—}• The matrix c ab is defined by 
c ++ = 1, c = —1, and c + ~ = c~ + = 0. The functional differentiation in Eq. ( |2.9| ) is 
carried out with variations in 5J + and 5J_ which satisfy the constraint that 5J + = 5J _ on 
the x° = x° hypersurface. The J± subscript in Eq. ( |2.9| ) indicates the functional dependence 
on J±, which has been shown to be causal )[28| . [27|| . In the limit J + = = J, the classical 

helds on the + and — time branches become equal, 



j+=j^=j 



j+=j-=j 


4>{x)j 


jJO, in|<f) H (x)|0, in)j, 


( 2 . 10 ) 


where 10, in) is the state which has evolved from the vacuum at to under the interaction 
'LhJ, and becomes the expectation value (^h) in the limit J = 0. The effective action is 
defined via the usual Legendre transform, with c ab now acting as a “metric” on the internal 
2x2 CTP held space, 


T[0+, 0_] = W[J+, J-} - c ab [ d 4 xJ a {x)Mx ), (2.11) 

J M 

where the J subscripts on ij>± are suppressed and the functional dependence of J± on 0 via 
inversion of Eq. (^.9| ) is understood. By direct computation, the inverse of Eq. ( |2.9| ) is found 
to be 


Ja{pC) 4>± Cab 


<*r[0+,0_] 

<50b (x) 


( 2 . 12 ) 


where we have indicated the explicit functional dependence of J± on 0± with a subscript, 
and c a b is the inverse of the matrix c ab defined above. In the limit 0 + = 0_ = 0, this yields 
the evolution equation for the expectation value j($ h )j = 0j in the state which has evolved 


3 Here, the choice of vacuum boundary conditions corresponds to adding a small imaginary part 
ie(0+ 2 — 0- 2 ) to the classical action 5 F . Alternatively, the boundary conditions correspond to 
the usual prescription m 2 —> m 2 — ie in 5' F [0], but with 5 ctp now redefined as 5 F [0 + ,0_] = 
S' F [0 + ] — S F [0_]*, where * denotes complex conjugation [221. 
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from |0, in) under the source interaction J<L H . The evolution equation for 0, the vacuum 
expectation value (0, in|$ H |0, in), is therefore 


hr[0 + ,0] 




hr[0 + ,0] 


0+— 4>-=4> 


5<f>_ 


= 0. 


(2.13) 


4>+— 0- =4> 


Using Eqs. (|2.12| ) and (|2.11| ), an integro-differential equation for T can be derived ||28 
in which the functional differentiations of T with respect to 0± are carried out with the 
constraint that the variations of (f>± satisfy 50+ = h0_ when x° = x®. The difference 0 a — 0 a 
is naturally interpreted as the fluctuations of a particular history 0 a about the “classical” 
held configuration (f) a . Let us, therefore, define the fluctuation field (p a = 0 a — 0 a or, in 
terms of Heisenberg held operators, 


Ph = $„ - (4„> = 4, 


(2.14) 


where angular brackets around the Heisenberg held operator <f> H denote an expectation 
value of $ H in the (time-independent) quantum state of the system. Performing the change 
of variables 0 a —> tp a in the functional integral, where 

[Pa = 0a 0a j (2.15) 


we obtain 


T[0+,0_] = -ih\n{ / Dif + Dip_ exp 

' ctp 


" i cf i 2 i i 5r[0 + ,0-] 

- S [0 + + ip + , 0— + ip-\ --y- <Pa 

n \ 00a 


This functional integro-differential equation has a formal solution |58 


A. A. j h A A. 

r[0+y_] = S p [<yy_] - —hidet^", 1 ) + r 1 [^ + .0_]. 


(2.16) 


(2.17) 


where x'), the second functional derivative of the classical action with respect to the 

held 0±, is 


h 2 S F 


iA ab (x, x') = 


S(j) a (x)(j) b (x 


A [0+104 


(2.18) 


The inverse of A ab is the one-loop propagator for the fluctuation held 0. The functional T | in 
Eq. ( j2.16| ) is dehned as — ih times the sum of all one-particle-irreducible vacuum-to-vacuum 
graphs with propagator given by A~^(x,x') and vertices given by a “shifted action” 5 F t , 
dehned by 


ip~] = 5 F b+ + 0 +, p- + 0 _] - s F [(/)+, 0 —] - / d 4 x 

Jm 

1 

~2 




S(pa 


d^x I d A x' 


' M 


IM 


5 2 S F 


54>a{x)54> b (x') 


■±] ) Pa 

77[0±] J Pa(x)cp b {x'). (2.19) 


7 

























For simplicity, we do not explicitly indicate the functional dependence of S? nt on <fi±. Fig. [I] 
shows the diagrammatic expansion for F]. Each vertex carries a spacetime label in M and 
a time branch label in {+,—}. The lowest-order contribution is order h 2 , i.e., at two loops. 
The propagator A~ x does not depend on h. The In det A term in Eq. (| 2 . 17 |) is the one- 
loop (order K) term in the CTP effective action. The CTP effective action, as a functional 
of 0±, can be computed to any desired order in the loop expansion using Eq. ( P~T71) . In 
general, this action contains divergences at each order in the loop expansion, which need to 
be renormalized. 

Functionally differentiating r[0 + ,0_] with respect to either 0 + or 0_ and making the 
identification 0+ = 0_ = 0 [as shown in Eq. ( |2.13| )1 yields a dynamical, real, and causal 
evolution equation for the mean held 0. Thus the 1PI effective action r[0±] yields mean- 
field, dynamics for the theory, which is a lowest-order truncation of the correlation hierarchy 
FlfTi] . However, for a detailed study of nonperturbative growth of quantum fluctuations 
relevant to nonequilibrium mean-held dynamics (or a symmetry-breaking phase transition), 
it is also necessary to obtain dynamical equations for the variance of <I> H , 


<* h > - (**r = (n) -r = (vi) = kg ++ (x,x\ 


( 2 . 20 ) 


where hG ++ (x,x') is the time-ordered Green function for the fluctuation held cp n , 
(T((p(x) n ip(x') H ))■ A higher-order truncation of the correlation hierarchy is needed in order 
to explicitly follow the growth of quantum huctuations; the 2PI effective action, to which 
we now turn, serves this purpose. 


III. TWO-PARTICLE-IRREDUCIBLE FORMALISM 


In a nonperturbative study of nonequilibrium held dynamics in the regime where quan¬ 
tum huctuations are significant, the 1PI effective action is inadequate because it does not 
permit a derivation of the evolution equations for the mean held ($ H ) and variance a t 

a consistent order in a nonperturbative expansion scheme. In addition, the initial data for 
the mean held 0 do not contain any information about the quantum state for huctuations 
tp around the mean held. The two-particlc-irreducible (2PI) effective action method can 
be used to obtain nonperturbative dynamical equations for both the mean held <p(x) and 
two-point function G(x,y), which contains the variance, as shown in Eq. ( |2.20|) . The 2PI 
method generalizes the 1PI effective action to an action T[0, G] which is a functional of 
possible histories for both 0 and G. Alternatively, the 2PI effective action can be viewed as 
a truncation of the master effective action to second order in the correlation hierarchy [03 


In this section we briefly review how the 2PI method works; more thorough presentations 
can be found in |32|,|3^|. 

Unlike the 1PI method where the mean held is fixed to be 0, the 2PI method Exes the 
mean held to be 0 and the sum of all self-energy diagrams to be G. This drastically reduces 
the number of independent diagrams which must be computed in order to obtain T[0,G] 
45| . Coupled dynamical equations for the evolution of 0 and G are obtained by separately 


varying T[0, G] with respect to G and 0. Imposing AT/A0 = 0 yields an evolution equation 
for the mean held 0, and setting 5Y/8G = 0 yields an evolution equation for G, the “gap” 
equation. The variance (<^|) is the coincidence limit of the two-point function hG, as seen 













from Eq. (|2.20|) . In a nonequilibrium setting, the closed-time-path method should be used 
in conjunction with the 2PI formalism in order to obtain real and causal dynamics for 0 and 
G |3|J3|,|451 . 

Let us apply the 2PI method to a scalar A0 4 theory in Minkowski space, with vacuum 
initial conditions. In a direct generalization of Sec. |I1]. both a local source J a (x ) and nonlocal 
source K a b(x,x') (which are c-number functions on M) are coupled to the held via hc ab J a (j) a 
and hc ab c cd K ac (x,x')(j)b(x)(j)d(x') interactions. Following Eq. 0 , the CTP generating func¬ 
tional is defined as a vacuum persistence amplitude in the presence of the sources J and K, 
which has the path integral representation 


Z[J, K\= / Dcj) + D(j)_ex p 

J ctp 


■^(s F [0 + ,0_] + [ d A xc ab J a (j)b 
n \ Jm 

+ 7 : [ d A x f dVc“ 6 c cd iP oc (a;,a; / ) 06 (a;) 0 d (a: / )') 
2 Jm Jm ' 


(3,1) 


Here, 5 F is as defined in Eq. ( j2.5|) , and we are using Z[J,K ] as a shorthand for 
Z[J + , J_; K++, K _, K + _,K_ + ]. The generating functional for normalized n-point func¬ 

tions (connected diagrams) is defined by 


The “classical” held ba(x) JK 

0 a (•£) jk 

fi>G a b(X) x ) JK 


W[J, K] = -ih\nZ[J,K}. 

and two-point function G a b(x,x') JK are then given by 
6W[J,K] 


C-ab " 


2CacCbd 


SJb(x) 

8W[J, K] 


6K cd (x,x') 


I \ ba(x^) JK (f)b(x ) j K t 


(3.2) 


(3.3a) 

(3.3b) 


where we use the subscript JK to indicate that 0 a and G a b are functionals of the sources 
J and K. 

In the limit K — J — 0, the classical held 0 a satishes 

(0+)j=X=0 = (0 -)j=A'=O = (0|$h|0) = 0; (3.4) 


i.e., it becomes the expectation value of the Heisenberg held operator <h H in the quantum 
state |0) (the mean held). In the same limit, the two-point function G a b is the CTP prop¬ 
agator for the fluctuation held defined by Eq. (|2.14| ). The four components of the CTP 
propagator are, for J = K = 0, 


hG ++ (x, x')\ J=K=0 = (0| T (tp H (x)p n (x')) |0), (3.5a) 

hG—(x,x')\ J=K=0 = ((f)\f (cp k(x)iPk(x')) 10), (3.5b) 

hG + -(x,x')\ J=K =o = (0 |¥ , h((c')^h((e)|0 ), (3.5c) 

hG_ + (x,x')\ J=K=0 = (0|<Ph(^)^h(^)I0)i (3.5d) 


in the Heisenberg picture. In the coincidence limit x' 
equivalent to the variance (<p|) defined in Eq. (j2.20p. 


= x, all four components above are 
Provided we can invert Eqs. (|3.3a|) 
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and (|3.3|) to obtain J and K in terms of 0 and G, the 2PI effective action can be defined as 
the double Legendre transform (in both J and K) of W[J,K] 


r [4>,G\ = W[J,K\- I d 4 xc ab J a (x)Mx ) 


' M 
1 


- / cl 4 x d 4 x'c ab c cd K ac (x,x')[hG bd (x,x') + 4> b (x)(p d (x')}. (3.6) 


IM 


i M 


As with W[J,K\, we are using r[0, G\ as a shorthand for r[0 + , 0_; G++, G _, G+_,G_ + ]. 

The JK subscripting of 0 and G is suppressed, but the functional dependence of 0 and G 
on J and K through inversion of Eqs. (|3.3a| ) and (|3.3| ) is understood. By direct functional 
differentiation of Eq. (b-6|), the inverses of Eqs.( |3.3a| ) and ( |3.3| ) are found to be 

^ ^ = -c ab J b (x) iG - \c ab c cd [ d 4 x'(K bd (x, x') iG + K db (x', x) iG )(f) c {x'), (3.7a) 

^ Jm 


S<j>a(x 

sr [ 0 ,G] 


h 


6G^(x,^ = -¥ CcUK ‘ AX ' X ' ha ’ 


(3.7b) 


where the subscript “0G” indicates that K and J are functionals of 0 and G. Once T[0,G] 
has been calculated, the evolution equations for 0 and G are given by 


ar[ 0 .G] 


H a (x) 

«h,G] 


= 0 , 


4>+— 4>~=4> 


5G ab (x,y) 


= 0 . 


(3.8a) 


(3.8b) 


Of course, the two equations contained in Eq. ( |3.8a|) are not independent, just as in 
Eq. (|2.13| ). In addition, only two of equations (|3.8b|) are independent, one on the diago¬ 
nal and one off diagonal in the “internal” CTP space. Using both Eq. ( |3.6| ) and Eq. ( |3.1|) , 
an equation for T[0, G\ in terms of the sources K and J can be derived, 


' ctp 


L h 


r[0, G\ = -ih\n{ I L)0 + L)0_exp 0 (5 F [0+, 0_] + c ab I d 4 xJ a (x)[(j) b (x) - (f> b (x)} 


IM 


,ac Jbd 


+-c ac d 


d x / d x'K ab (x,x')[(j) c (x)(j) d (x') - cp c (x)(j) d (x') - hG cd (x,x')\ 


' M 


' M 


(3.9) 


The sources K and J in the right-hand side of Eq. ( |3.9| ) are functionals of 0, through 
Eqs. ( p.7a| ) and (|3.7b| ). Expressing this functional dependence, we obtain a functional 
integro-differential equation for T, 


r [ 0 ,G] = d 4 x d 4 x 


4 , W[0,G] 


Jm 
— ill In 


M 6G ba (x',x ) 


Gab(x,x') 


’ ctp 


D<j) + D(f)_ exp ^(S F [0 + ,0_] - / "~' J (0 O - 0 O ) 


TT[ 0 ,G] 


i 


d 4 x / d 4 x' 


4 , 5r[0,G] 


' M 


M 5G ba (x’,x) 


JM 50 a 

[ 0 a (%) ~ </> a ( x )][M X ') “ M X ')\ 


(3.10) 
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We have dropped the JK subscripting because the functional derivatives in the equation are 
only with respect to 0 and G. As in Sec. |J, a change of variables D(j)± —»• Dip± is carried 
out in the functional integral. The resulting equation 


,4, «y,G] 


r IM= d'G } 

JM JM V(j ba \X i -L) 


- rG a b(x,x') 


-ih\ n<j j D(p + Dip_ex p - ( 5 F [<^+ + <f> + , ip- + 0_] 


^4 ^ '-'J J- / / J4„/ ^r[0, G] 


'M 


fEhihl. _ I [ d * x I dV 

50 a ftJM Jm 5G ba (x',x ) 


f (3-11) 


has the formal solution 32 


„ a, jh 

T[0, G] = £ F [0] ——lndet(G a (,) 

ih 


+ — / d 4 x / d A x A ab (x ) x)G ab (x, x) + r 2 [0, G], (3.12) 

4 Jm Jm 


where _4. a6 is the second functional derivative of the classical action G F , evaluated at 0 a . 
The functional T 2 is —ih times the sum of all two-particlc-irreducible vacuum-to-vacuum 
diagrams with lines given by G ab and vertices given by a shifted action £ F t . We have, for 


A 


ab 


and, for S F 

' in 


<5 2 £ F 


iA ab (x, x') = 


8J) a {x)J) b {x') 


^intb] 



1 

2 



' M 


d A x'A ab (x, x')ip a (x)(p b (x'). 


The shifted action for the A0 4 scalar held theory is 


(3.13) 


(3.14) 


$L b] = S L b+] - s lt b-] , ( 3 - 15 ) 

in terms of 

SfM = J M d * x (j<b + , (3-16) 

where the functional dependence of G F nt on <j>± is not shown explicitly. Two types of vertices 
appear in Eq. (|3.16|) : a vertex which terminates four lines and a vertex terminating three 
lines which is proportional to the mean held 0. The expansion for T 2 in terms of G and 
0 is depicted graphically up to three-loop order in Fig. |2j. Each vertex carries a spacetime 
label in M and time branch label in {+, —}. In general, the 2PI effective action contains 
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divergences at each order in the loop expansion. It has been shown formally that if the held 
theory is renormalizable in the “in-out” formulation, then the “in-in” equations of motion 

IB1. 


are renormalizable 


In the closed-time-path formalism it is easier to carry out explicit 


renormalization in the equations of motion, i.e., the mean-held and gap equations, which we 
will do in Sec. IVT5L 


Various approximations to the full quantum dynamics can be obtained by truncating the 
diagrammatic expansion for IV Throwing away ^ in its entirety would yield the one-loop 
approximation. In Fig. |2|. there are two two-loop diagrams, the “double bubble” and the 
“setting sun.” Retaining just the double-bubble diagram yields the time-dependent Hartree- 
Fock approximation fl32j. Retaining both diagrams gives a two-loop approximation to the 


theory.^ This approximation will yield a non-time-reversal-invariant mean-held equation 
above threshold, due to the setting sun diagram |44|] . The time-reversal noninvariance of the 
mean-held equation generated by the 2PI effective action is a consequence of the fact that 
the 2PI effective action really corresponds to a further approximation from the two-loop 
truncation (in the sense of topology of vacuum graphs) of the master effective action |44| . 
The two-loop truncation of the master effective action is a functional Ti = 2 [(f>, G, C 3 ] which 
depends on the mean held 0, the two-point function G, and the three-point function C 3 . 
The four-point function C 4 also appears, but is not dynamical due to a constraint. The full 
set of equations, 


5T l= S,G,C 3 ] 

Ha 

5T l=2 [4>,G,C 3 } 

5G ab 

hr«= 2 [ 0 , g, c 3 \ 

5(C 3 )abc 


0, 

0, 

0, 


(3.17a) 

(3.17b) 

(3.17c) 


is time-reversal invariant. However, the 2PI effective action is obtained by solving Eq. ( |3.17c|) 
with a given choice of causal boundary conditions and substituting the resulting C 3 into T; = 2 , 
to obtain T 2 [0, G]. This “slaving” of C 3 to 0 and G with a particular choice of boundary 
conditions is what breaks the time-reversal invariance of the theory |F|]. In the paper fol¬ 
lowing |5(J where we discuss the preheating dynamics, we work with further approximations 
which discard the setting sun diagram, and thus regain time-reversal-invariant equations. 


IV. A $ 4 FIELD THEORY IN CURVED SPACETIME 

In this section the quantum dynamics of a scalar A $ 4 held theory is formulated in semi- 
classical gravity, where the matter fields are quantized on a curved classical background 


4 A different approximation, the 1/A expansion, is used in Sec. |V] to study the nonequilibrium 
dynamics of the 0(A) field theory. 
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spacetime.[] The two-particle-irreducible effective action is used in conjunction with the 
CTP formalism to obtain coupled evolution equations for the mean field (<h H ) and variance 
(<f)|) — (<f> H ) 2 in the A<F 4 model which are manifestly covariant. 

Let us consider a quartically self-interacting scalar field 0 in a globally hyperbolic, curved 
background spacetime with metric tensor g^ u . The diffeomorphism-invariant classical action 
for this system is 


s[h<n = s a \ ! n+&[<!>, gr\. 


(4.1) 


where g^ u is the contravariant metric tensor, and S G and S F are the classical actions of the 
gravity and scalar field sectors of the theory, respectively. For the scalar field action, we 
have 


S F [(ft,g 


flU] 



+ m 2 + (R)& + 



(4.2) 


where £ is the (dimensionless) coupling constant to gravity (necessary in order for the field 
theory to be renormalizable |)9|), □ is the Laplace-Bcltrami operator in terms of the covari¬ 
ant derivative V M , and R is the scalar curvature. The constant m has units of inverse length, 
and the (ft self-coupling A has units of 1/h. Following standard procedure in semiclassical 
gravity ||2|, we define the semiclassical action for gravity to be 

sG [gn = Jd^Xyf^g [R- 2A + cR 2 + bR a0 R aP + aR a ^ 5 R a ^\ , (4.3) 


where a, b, and c are constants with dimensions of length squared, Rafiks is the Riemann 
tensor, R a p is the Ricci tensor, A is the “cosmological constant” (with units of inverse length- 
squared), y/^g is the square root of the determinant of g )iu , and G (with units of length 
divided by mass) is Newton’s constant. As a result of to the generalized Gauss-Bonnet 
theorem |60|, the constants a, b, and c are not all independent in four spacetime dimensions; 
let us, therefore, set a = 0. Classical Einstein gravity is obtained by setting 6 = 0 and c = 0. 
Minimal and conformal coupling (for the (ft field to gravity) correspond to setting £ = 0 and 
£ = 1/6, respectively. 

The motivation for including the arbitrary coupling £ and the higher-order curvature 
terms R 2 and R al3 R a /3 in the classical action S is that we wish to study the semiclassical 
dynamics of the theory. In the semiclassical gravity field equation and matter held equations, 
divergences arise which require a renormalization of b, c, G, A, m, £, and A pi]. These 
quantities are understood to be bare; their observable counterparts are renormalized. 

The classical Euler-Lagrange equation for (ft is obtained by functionally differentiating 
S F [(ft,g flu ] with respect to (ft , and setting 5S F /S(ft = 0, 


5 The semiclassical approximation is consistent with a truncation of the quantum effective action 
for matter fields and gravity perturbations at one loop [i.e., at order 0{R)\ |29|] or [in the case of 
the O(N) field theory studied here] at leading order in the 1/N expansion ^5 [|r]]. 
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□ + m 2 + 


( 4 . 4 ) 


The Euler-Lagrange equation for the metric g^ y is obtained by functional differentiation of 
S with respect to (it is assumed that the variations S(p and Sg lliy are restricted so that 
no boundary terms arise), 

Gfiu + Ag^ v + c + b ® H^ u = SttGT^, (4-5) 


where the tensors G^, 
energy-momentum tensor, 


and are 


defined by |UT|JU^] and TJ 


[MS 


is the classical 


T,Pv = (1 - + ^ 2 ^ - 0 g»v9 pc <p-p(p-a - 2 + 2£g ^Ucj) 


-iG^cp 2 + ^g, w Cm 2 + ^(p 2 ' ^ 


(4.6) 


We are interested in the dynamics of expectation values in the semiclassical theory, which 
in nonequilibrium field theory does not follow directly from functional differentiation of 
the usual Schwinger-DeWitt or “in-out” effective action. Instead, the Schwinger-Keldysh 
formalism (reviewed in Sec. p) should be used. Here we discuss the implementation of the 
Schwinger-Keldysh method in curved spacetime. 

The first step is to generalize the closed-time-path (CTP) manifold A4, defined in 
Eq. (j2.3|), to curved spacetime. Let E* be a Cauchy hypersurface chosen so that its past 
domain of dependence |k|, _D_(E*), contains all of the dynamics we wish to study. Let us 


then define the manifold (with boundary) 


M = £>_(£*)■ (4.7) 

The CTP manifold Ai is defined following Eq. Q2.3|) as a quotient space constructed by 
identification on the hypersurface E* C DM as in Eq. ( |2.3| ) where the equivalence relation 
is the same as Eq. ( |2.4|) except that the matching of + and — time branches is now done on 
E*. We construct an orientation on Ai using the canonical volume form from M, Em, and 
define the volume form on Ai to be 


f e M on M x {+}, 
I - e m on M x {-}. 


(4.8) 


Finally, we let (p and g^ v be independent on the + and — branches of A4, provided that 
g+ u = g'A and <p + = (p _ on E*. In other words, (p and g tw must be a scalar and a tensor, 
respectively, on Ai. In terms of the volume form Em, we can write a scalar field action on 

M, 


S F [<P±,g*n = S F [</>+, g'p 1 '] - S*[<P-,g?l, (4-9) 

where S F [(p] is given by Eq. (|4.2|) , and g± is the metric tensor on the + and — branches of 
Ai. Using Eq. (|4.3|) we can similarly define the gravity action S G on Ai, 
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s a \tf,<F] = s a \g!?]-s a \ ! n, (4-io) 

where it is understood that only configurations of g± satisfying the constraint gffi = gf u on 
E* are permitted. 

In semiclassical gravity the scalar field theory (with action S F ) is quantized on a classical 
background spacetime, with metric g^, whose dynamics is determined self-consistently by 
the semiclassical geometrodynamical field equation. Let us denote the Heisenberg-picture 
field operator for the canonically quantized <f field by <f> H . We wish to compute the quantum 
effective action F for this scalar field theory, using the two-particlc-irreducible (2PI) method 
described in Sec. B In terms of S F (now defined on the curved manifold A!), we define 
the 2PI, CTP generating functional Z[J, K , g^] as follows: 


Z[J,K,g^} = / D(f) + D(f)_ exp 


' ctp 


h 


(S P [£fc,0f] + [ d 4 Xy/=f c C abc J c 
V J M 


+- I d A xy/-g a i j d A x \f—g l c ,c aba ' d :dc ' K ac (x,x')(i> b (x)<p d (x') 


IM 


'M 


(4.11) 


where we have written Z[J, K, g^] as a shorthand for Z[J±, K±±, g±]. The three-index 
symbol c abc is defined by 


c abc 


1 if a = b = c = +, 
—1 if a = b = c = —, 
0 otherwise. 


(4.12) 


The boundary conditions on the functional integral define the initial quantum state (as¬ 
sumed here to be pure). In this and a subsequent paper (in which preheating dynamics 
of inflationary cosmology is studied ||5D[), we are interested in the case of a quantum state 
corresponding to a nonzero mean field 0, with vacuum fluctuations around the mean field. 
This entails a definition of the vacuum state for the fluctuation field, defined in Eq. (|2.14|) . In 
curved spacetime in general, there does not exist a unique Poincare-invariant vacuum state 
for a quantum field [0]||. For an asymptotically free field theory, a choice of “in” vacuum 
state corresponds to a choice of a particular orthonormal basis of solutions of the covariant 
Klcin-Gordon equation with which to canonically quantize the field. 

From Eq. ( |4.11| ), we can derive the two-particlc-irreducible (2PI) effective action 
r[0, G, g^} following the method of Sec. |TJ. with the understanding that F now depends 
functionally on the metric g± on the + and — time branches. The functional differentiations 
should be carried out using a covariant generalization of the Dirac 6 function to the manifold 
M |2|]. The functional integro-differential equation ( |3. 1 ID for the CTP-2PI effective action 
can then be generalized to the curved spacetime M in a straightforward fashion, modulo the 
curved-spacetime ambiguities in the boundary conditions of the functional integral (|4.11|) . 

The (bare) semiclassical field equations for the variance, mean field, and metric can then 
be expressed in terms of variations of S G [g lw ] + F[fi, G,g^ v ] with respect to G±±, <f>±, and 
gff , respectively, followed by metric and mean-field identifications between the + and — 
time branches, 


VM + riic.ir]) 




fiv 

a 


= o, 


(4.13a) 


7 7 7 U V UV ,,,, 

<(> + =</>_= 0 ; gy =gt ; =g^ 
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ST[lG,g^\ 


5 (fia 

ST[lG,g^] 


5G, 


ab 


7 7 7 U.V UV . n , 

<f>+ = <P- =0; 5+ =9- =9^ 


777 /XI/ ill/ 

0+=<A_=0; 3+ =31 =3^ 


= 0 , 


= 0. 


(4.13b) 

(4.13c) 


As above, CTP indices are suppressed inside functional arguments. Eqs. (|4.13| a), (|4.13| b), 
and ( |4.13| c) constitute the semiclassical approximation to the full quantum dynamics for the 
system described by the classical action S' G [g / “'] + S' F [</>, g^}. The semiclassical held equation 
(with bare parameters) for g ,JV is obtained directly from Eq. (|4.13|a), 


+ A gpv + c + b = -8nG(T^), 

where (T /IU ) is the (unrenormalized) energy-momentum tensor defined by 


(r„„) = 




V -9 v % 




(4.14) 


(4.15) 


*+=,£_=*; =y_ 


r=a^=a^ 


Equation ( 4.14 ) gives the spacetime dynamics; the dynamics of cp and G are given by the 
mean-held and gap equations, (|4.13| b,c). In Eq. (|4.15|) , the angle brackets denote an ex¬ 
pectation value of the energy-momentum tensor (with the Heisenberg held operator <h H 
substituted for cp in the classical theory) with respect to a quantum state \(p) defined by the 
boundary conditions on the functional integral in Eq. ( |4.11| ). In four spacetime dimensions 
(unrenormalized) (T /lu ) has divergences which can be absorbed by the renormalization of G, 
A, b, and c [j2]. This renormalization should be carried out in the held equations rather than 
in the CTP effective action ||28|| . 

The energy-momentum tensor as defined in Eq. ( |4.15|) is obtained by variation of the 
2PI effective action T, which is a functional of the metric g± on both the + and — time 
branches. From Eq. (|4.11|) , it is possible to derive T[(p,G, g^] as an arbitrary functional of 
g^ and g^. However, in practice it is often easier to work in the simplified case where the 
metric is hxed to be the same on both the + and — time branches, i.e., 


g+ = 9- = 9^, 


(4.16) 


in the computation of T[0, G, g ,JiV ]. Once T[0, G, g ,lv ] (or some consistent truncation of the full 
quantum ehective action for 5 F ) has been computed, it is then straightforward to determine 
how T[0, G, g^] should be generalized to the case of an arbitrary metric on Ad, for which g+ u 
and g*py are independent. The bare energy-momentum tensor (T^ u ) can then be computed 
using Eq. ( |4.15| ). Accordingly, in Sec. 0, we fix g'+ = g^ = g^ u in the calculation of 

ro.G.g""}. 

The semiclassical Einstein equation is a subcase of the general geometrodynamical held 
equation ([4.14 ), obtained (after renormalization) by setting the renormalized b = c = A = 0 
(assuming no cosmological constant) 0: 


G^ = - 8ttG(T p v ). 


(4.17) 
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Having shown how to derive coupled evolution equations for the mean field, variance, and 
metric tensor in semiclassical gravity, we now turn our attention to the scalar O(N) model 
in curved spacetime. 


V. O(N) FIELD THEORY IN CURVED SPACETIME 


In this section we derive coupled nonperturbative dynamical equations for the mean held 
0 and variance for the minimally coupled O (N) scalar held theory with quartic self¬ 
interaction and unbroken symmetry. The background spacetime dynamics is given by the 
semiclassical Einstein equation. These equations take into consideration the back reaction 
of quantum particle production on the mean held, and quantum helds on the dynamical 
spacetime, self-consistently. In our model the Heisenberg-picture quantum state |0) is a 
coherent state for the held $ H at the initial time r] 0 , in which the expectation value (<f> H ) is 
spatially homogeneous. The coherent state is defined with respect to the adiabatic vacuum 
constructed via matching of WKB and exact mode functions for the fluctuation held in some 
asymptotic region of spacetime. 

The O(iV) held theory has the property that a systematic expansion in powers of 1 /N 
yields a nonperturbative reorganization of the diagrammatic hierarchy which preserves the 

M 


Ward identities order by order 


Unlike perturbation theory in the coupling constant, 


which is an expansion of the theory around the vacuum configuration, the 1/N expansion 
entails an enhancement of the mean held by \/N ; this corresponds to the opposite limit 
of strong mean held. (This is precisely the situation which can arise in chaotic inflation 
at the end of the slow-roll period, where the inhaton mean held amplitude can be as large 
as Mp/3 ||50|.) As discussed in Secs. |Hj and ES the nonequilibrium initial conditions for 
the mean held as well as the nonperturbative aspect of the dynamics requires use of both 
closed-time-path and two-particle-irreducible methods. The 1/N expansion can be achieved 
as a further approximation from the two-loop, two-particle-irreducible truncation of the 
Schwinger-Dyson equations. 

Although in this study we assume a pure state, the 2PI formalism is also useful for an 
open system calculation, in which the mean held is defined as the trace of the product of 
the reduced density matrix p and the Heisenberg held operator $ H , Tr(p$ H ). When the 
position-basis matrix element (</>i|p(?7o)|02) can be expressed as a Gaussian functional of 
and 02j the nonlocal source K can encompass the initial conditions coming from p(t 0 ) 

m 


in a natural way 


In order to incorporate a density matrix whose initial condition is 


beyond Gaussian order in the position basis, one can work with a higher-order truncation 
of the master effective action & The leading-order 1/N approximation is equivalent to 
assuming a Gaussian initial density matrix; therefore, the 2PI effective action is adequate 
for our purposes. 


A. Classical action for the O(iV) theory 

The O (N) held theory consists of N spinless helds 0 = {0*}, i = 1,..., N, with an action 
which is invariant under the AT-dimensional real orthogonal group. The generally covariant 
classical action for the O(iV) theory (with quartic self-interaction) plus gravity is given by 
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S[#,9'‘1 = S g [9'‘1 + S'[^,9' , 1, 


( 5 . 1 ) 


where S G [g^ u ] is defined in Eq. ([4.3| ) for the spacetime manifold M with metric g^, and the 
matter field action S F [0®, g^} is given by 


S F [(j)\g^} = ~2 I dix V~9 


' M 


$-(D + m 2 +(R)$+^-f ) 2 


(5.2) 


The 0(1V) inner product is defined by[] 


(£■(£= 4> l (jy i 8i j . 


(5.3) 


In Eq. (|5.2|) , A is a (bare) coupling constant with dimensions of 1 /h, and £ is the (bare) 
dimensionless coupling to gravity. The classical Euler-Lagrange equations are obtained by 
variation of the action S separately with respect to the metric tensor g^ and the matter 
fields q'j l . In the classical action (|5.2|) , the O(IV) symmetry is unbroken. However, the O(IV) 
symmetry can be spontaneously broken, for example, by changing m 2 to — m 2 in S F . In the 
symmetry-breaking case with tachyonic mass, the stable equilibrium configuration is found 
to be 


2Nm o 

• 0 = —= V 2 , 


(5.4) 


which is a constant. If we wish to study the action for oscillations about the symmetry- 
broken equilibrium configuration, the O(IV) invariance of Eq. (|5.2|) implies that we can 
choose the minimum to be in any direction; we choose it to be in the first, i.e., (0 1 ) 2 = v 2 . 
In terms of the shifted field a = 0 1 — v and the unshifted fields (the “pions”) 7r* = q i*, 
z = l,..., — 1, the action becomes 


S F [(T,Tr,g liv } = -~ / d^x^g 


IM 


(x(D + m 2 + iR)a + 7? • (□ + m 2 + £R)n + 2(m 2 + £,R)a 2 

(5.5) 


+2 \ —Ma 3 + 2\ — Mtt ■ ncr H — cr 4 - if • ncr 2 H — (if • if) 2 

V 2 V 2 4 2 4 V 


One can show that the effective mass of each of the “pions” if (defined as the second derivative 
of the potential) is zero, due to Goldstone’s theorem. The theorem holds for the quantum- 
corrected effective potential as well [Q. In this paper we study the unbroken symmetry 
case, in order to focus on the parametric amplification of quantum fluctuations; this avoids 
the additional complications which arise in spontaneous symmetry breaking, e.g., infrared 
divergences ||T9| , |65| , |55i . 


6 In our index notation, the latin letters i,j,k,l,m,n are used to designate 0(N ) indices (with 
index set {1,..., N}), while the latin letters a, b, c, d, e, f are used below to designate CTP indices 
(with index set {+,—}). 
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B. Quantum generating functional 


We aim at deriving the mean-field and gap equations at two-loop order. The 2PI generat¬ 
ing functional for the 0 (N) theory in curved spacetime is defined using the closed-time-path 
method in terms of c-number sources J* and nonlocal c-number sources K^ b on the CTP 
manifold A4. 




i,a 


' ctp 


h 


T\S F [^,g^}+ / d A x^gc ab J c 


' M 


+ - / d A x^g / d 4 x'^g'c ab c cd K l a J c (x,x')(l)^x)(j) l d (x , )S ik S j i 

^ -/M JM 


, (5.6) 


where the CTP classical action is defined as in Eq. ( [i.9p , with 0! t replacing <J)±, and the 
time branch indices on g fW suppressed. The sources J* are coupled to the field by the O (N) 
vector inner product 


Jn 


= w 


V 


(5.7) 


The time branch labels on the metric tensor are suppressed for simplicity of notation. [| 
The designation “CTP” on the functional integral indicates that we sum only over field 
configurations for (j) l a on M. which satisfy the condition 0+ = 0!_ on £*, where £* is defined 
in Sec. E3 In addition, the boundary conditions on the asymptotic past field configurations 
for 0! t in the functional integral correspond to a choice of “in” quantum state |0) for the 
system. The generating functional for normalized expectation values is given by Eq. ( |3.2|) , 
with the additional functional dependence of both W and Z on understood. As above, 
we henceforth omit all indices in functional arguments. In terms of this functional, we can 
define the “classical” field 0 and two-point function G by functional differentiation, 



^ 8W sa 

V=gdJl(x) ’ 


(5.8) 


fta( X )<K( X ') + hG ab( X : X ') = 2 ‘ 


Cbd 


SW 


V~9 \/w7 $K l c %(x,x') 
In the zero-source limit K^ b = J l a = 0, the classical field (p l a satisfies 

(0V)j=iC=o = (0 ? _)j=A'=0 = (0|$hI0) = 0* 


5 ik 8 jl . 


(5.9) 


(5.10) 


‘The suppression of CTP indices on the metric tensor does not prevent computation of the ex¬ 
pectation value of the stress-energy-momentum tensor. It will be clear how to reinstate the time 
branch indices after the 2PI effective action has been explicitly computed, e.g., in the large-A 
approximation. 
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as an expectation value of the Heisenberg field operator in the quantum state |0). The 
fluctuation field is defined [as in Eq. (|2.14|) 1 in terms of the Heisenberg field operator <h H and 
the mean field 0 (times the identity operator), 

■rf, = *•„-#■ (5.11) 

In the same limit J — K — 0, the two-point function G^ b becomes the CTP propagator for 
the fluctuation field. The four components of the CTP propagator are (for J* = K^ b = 0) 


hG+ + ( x , x') \j =K= o = (0|T((p[ I (a;)(p4(a; / ))|0), 

(5.12a) 

hG^_(x,x')\ J=K=0 = (0|T(<p® H (a;)<^(a/))|0), 

(5.12b) 

hG l l_(x,x ') J=K =o = (0|^(x / X(:r)|0), 

(5.12c) 

hG l l + (x,x')j =K =o = (0|/ H ((r)^((r / )|0). 

(5.12d) 


In the coincidence limit x' = x, all four components above are equivalent to the mean-squared 
fluctuations (variance) about the mean field 0®, 

hGl + {x,x)\ J=K =o = (0|(<4) 2 |0> = (Oh) 2 )- (5-13) 

Provided that the above equations can be inverted to give J® and K^ b in terms of <f> l a and 
G^b, we can define the 2PI effective action as a double Legendre transform of W, 


T[t>,G,r} = W{J,K,r]- / d‘xV=9^J'M 


IM 


*7 


1 

2 


'M 


d^x^g / hG b l d (x,x') + 4> b (x)4> l d (x') 5 ik 5 jh (5.14) 


' M 


where J* and K^ b above denote the inverses of Eqs. (|5.8|) and (|5.9| ). From this equation, it is 
clear that the inverses of Eqs. (|5.8|) and (|5.9|) can be obtained by straightforward functional 
differentiation of T, 

= c ab Sij (-0 0) - ^(lVv C ?[4('h x> ) + K db( x ', x )]<t> l dS k i) ( 5 - 15 ) 


V-9S(f>i(x) 


1)1 1 - 7 = -^ caCcbd Kt d (x',x). 


V=9dGZ(x,x')V=g 


(5.16) 


Performing the usual field shifting involved in the background field approach |5S|, it can 
be shown that the 2PI effective action which satisfies Eqs. (|5.14j) , ( |5.15|) , and (|5.16|) can be 
written 

T\i,G,r\=S r [i>,r] - ylndet [G«] 

+ V l d>x 'f~s [ 1)6'“;!.,- r 2 y, G.if"'', (5.17) 

* J M J M 
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where the kernel A is the second functional derivative of the classical action with respect to 
the held 0, 


iAij(x,x') = 


( 


8 2 S F 




77 y>,9- 


(IV 


v=?’ 


(5.18) 


and T 2 is a functional to be defined below. Evaluating Af- by differentiation of Eq. 
we hnd 


iAij(x, x')= Si jC ab [O x + m 2 + £R(x)] 


n abcd 


2N 


[4 > c( x )$' d (x)}5 ij 8 k i + 24> k c (x)ft d (x)5 ik 8ji )<5 4 (x - x') 




(5.19) 


where the four-index symbol c abcd is dehnecl in exact analogy with Eq. ( |4.12|) . In Eq. ( |5. 1 7D , 
To is — ih times the sum of all two-particle-irreducible vacuum-to-vacuum graphs with prop¬ 
agator G and vertices given by the shifted action «S? t , defined by 




— - / d A x I d^x' 


IM JM 


8 2 S F 


d 4 x 


IM 


8S F 

Scpi 


A 


flU 


Va 




(5.20) 


The expansion of T 2 in terms of G and 0 is depicted graphically in Fig. |2|. From Eqs. (|5.20|) 
and (U), S F t is easily evaluated, and we hnd 


5■>, <r] = sl [<e + ,<n - s, F >-, <H, 

in terms of an action S F nt on M dehned by 


sr>,<n = 


A 


2N 


d^xsj—g 


'M 


-(0-0) 2 + (4> ■$)($■$) 


(5.21) 


(5.22) 


The two types of vertices in Fig. |2| are readily apparent in Eq. (|5.22|) . The hrst term 
corresponds to the vertex which terminates four lines; the second term corresponds to the 
vertex which terminates three lines and is proportional to 0. 

The action T including the full diagrammatic series for T 2 gives the full dynamics for 
0 and G in the O (N) theory. It is of course not feasible to obtain an exact, closed-form 
expression for T 2 in this model. Various approximations to the full 2PI effective action can 
be obtained by truncating the diagrammatic expansion for T 2 . Which approximation is most 
appropriate depends on the physical problem under consideration. 

1. Retaining both the “setting-sun” and the “double-bubble” diagrams of Fig. || corre¬ 
sponds to the two-loop, two-particle-irreducible approximation [44|. This approxima¬ 
tion contains two-particle scattering through the setting-sun diagram. 

2. A truncation of T 2 retaining only the “double-bubble” diagram of Fig. |] yields equa¬ 

tions for 0 and G which correspond to the time-dependent Hartree-Fock approxima¬ 
tion to the full quantum dynamics This approximation does not preserve 

Goldstone’s theorem, but is energy conserving (in Minkowski space) 
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3. Retaining only the (tr G^ b ) 2 piece of the double-bubble diagram corresponds to taking 
the leading order 1/iV approximation, shown below in Sec. |V Q 

4. A much simpler approximation consists of discarding lb altogether. This yields the 
one-loop approximation, whose limitations have been extensively documented in the 
literature fTT].|To].[TT].[T^| . 

Let us first evaluate the 2PI effective action at two loops . This is the most general 

of the various approximations described above. Here both two-loop diagrams in Fig. are 
retained. The 2PI effective action is given by Eq. (|5.17|) , and in this approximation, T 2 is 
given by 

A h 2 


Ti[‘t>,G, s r\ = 


4JV 


—c 
2 


,abcd 


d 4 x V=g{GZ(x, x)G%(x, x) + 2 G%(x, x)G J J d (x, x)]SijS kl 


IM 


iX 

'N 


(flbcd ^fl! b' c! d! 


IM 


d A x\f^g / d A x'^^^{x)^(x')[Gw(x : x')G J c {,(x,x')G k d d,(x,x') 


IM 


+2 G^,(x, x')G 3 £(x, x')G%(x, x')] 8 jk 5 jlk , 


(5.23) 


Functional differentiation of Y with respect to 0 and G leads to the mean-held and gap 
equations, respectively. The gap equation obtained at two loops is given by 

(G- 1 )£(x,x') = + l ^c M S\x - x') [6‘iS M G“(x,x) + 2 G» d (x,x)\ 


h\ 2 

2N 2 


^acde^bc'd'e' 


dkk’div 


k(x)^A x, )GM'{x, x') G k J(x, x') 

+20c ( x )(j) l c i(x')G dd ! (x, x')G l e J e ,{x, x') 
+2 ^ c (x)^(x')G^ d ,(x, x')G 2 '(x, x') 

+2 |(x)^(0)Gav')G^(x,0) 

+2 4> k c (x)& c ,(x')Gfy(x, x')G 3 £(x, x') 


The mean-held equation is found to be 

+ rn 2 + (R) + 


ih 2 x 2 


E cm (x,x') 


4 N 2 JM 

fc \„abcd , 

~ 2 ir{ 8 ijW<%(x, x) + S^ l d [GZ(x, x) + G^(x, x 


in terms of a nonlocal function T, cm (x, x’) defined by 

„/\ ~ebcd jj! b'c!d! 2i ( __/A s~imi' / „ ™ 


Y, em (x, x') = c ebcd c a 


G$ (x,x')G%(x,x')G%(x,x') 

+ 2 G T/{ X , x ') gJ cc>( x , n')G$' d (x, x') 
+G^{x\x)G^x\x)Gf d {x\x) 

+ 2 GZ(x', x)G k X(x ', x)G^ b (x', x) 5 jk 5 jlk ,5 iV . 


(5.24) 


(5.25) 


(5.26) 
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Taking the limit <p l + = 0L = 0* in Eqs. (|5.24|) and (|5.25|) yields coupled equations for the 
mean held 0* and the CTP propagators G^ b , on the fixed background spacetime g^. The 
equations, as well as the semiclassical Einstein equation obtained from Eq. ( 4.13| a), are real 
and causal, and correspond to expectation values in the 0 + = 0 _ = 0 limit. The 0( A 2 ) parts 
of the above equations are nonlocal and dissipative. The nonlocal aspect makes numerical 
solution difficult; the dissipative aspect will be addressed in a future publication |66|]. One 


can regain the perturbative (amplitude) expansion for the CTP effective action at two loops 
by expanding the one-loop CTP propagators in Eq. ( 5.25 ) in a functional power series in 0 . 


C. Large-IV approximation 


We now carry out the 1/iV expansion to obtain local, covariant, nonperturbative mean- 
held and gap equations for the O(iV) held theory in a general curved spacetime. The l/N 
expansion is a controlled nonperturbative approximation scheme which can be used to study 
nonequilibrium quantum held dynamics in the regime of strong quasiclassical held amplitude 
In the large- N approach, the large-amplitude quasiclassical held is modeled by 
N helds, and the quantum-held-theoretic generating functional is expanded in powers of l/N. 
In this sense the method is a controlled expansion in a small parameter. Unlike perturbation 
theory in the coupling constant, which corresponds to an expansion of the theory around 
the vacuum, the large-iV approximation corresponds to an expansion of the held theory 
about a strong quasiclassical held configuration J3(J. At a particular order in the l/N 
expansion, the approximation yields truncated Schwinger-Dyson equations which are gauge 
and renormalization-group invariant, unitary, and (in Minkowski space) energy conserving 
|36| . In contrast, the Hartree-Fock approximation cannot be systematically improved beyond 
leading order and (in the case of spontaneous symmetry breaking) violates Goldstone’s 
theorem [f4t|. 

Let us implement the leading order large-iV approximation in the two-loop, 2PI mean- 
held and gap equations (|5.25|) and (|5.24|) derived above. This amounts to computing the 
leading-order part of T in the limit of large N, which is O(N). In the unbroken symmetry 
case, this is easily carried out by scaling 0 by y/N and leaving G unsealed ||32|1 , 


4>a( x ) 

Gab(x,x') 

ti(x) 


VN(f)a(x), 

G ab (x,x')8 ij , 

Vaix), 


(5.27a) 

(5.27b) 

(5.27c) 

(5.27d) 


for all i,j. The Heisenberg held operator (p l H scales like (p l a in Eq. (|5.27d|) . In the above 
equations, the connection between the large -N limit and the strong mean-held limit is clear. 

The truncation of the l/N expansion should be carried out in the 2PI effective action, 
where it can be shown that the three-loop and higher-order diagrams do not contribute (at 
leading order in the l/N expansion). Let us now also allow the metric g^ u to be specihed 
independently on the + and — time branches. We find, for the classical action, 


S^g^} = S^g^}- S^g^ 


(5.28) 


23 





























where 


S F \(j),9 fl 1 = ~ [ d^x^g 

1 J M 

The inverse of the one-loop propagator is^ 


0(D + m 2 + £R)0 + -0 


iA ab (x,x') = - 


X 


c abc [ □* + m 2 + £,R c (x)] + -c abcd <j> c (x)Mx) 


5 4 (x — x') 


\f-g i b 


(5.29) 


(5.30) 


Finally, for the CTP-2PI effective action at leading order in the 1/N expansion, we obtain 

jhN 

r= S r [4>,g» 1 - —Indet [G afe ] 


iftlV 

~ J 

Xh 2 N 


d A x^g a / dV v / ^ 7 fe-4 ab (ic / ,ic)G afe (a:,a; / ) 


M 

_i abode 


IM 


d*xy/^g e G ab (x, x)G cd (x, x) + 0(1). 


(5.31) 


<M 


Applying Eq. ( |4.13| b) and taking the limits 0 + = 0_ = 0 and g'+ = g'A = g we obtain 
the gap equation for G ab at leading order in the 1/N expansion, 


{G~ l ) ab {x,x') = A ab (x,x’) + l ^c abcd G cd {x, x)8\x - x')-jL= + O , (5. 


32) 


where 


iA ab (x, x') = — 


c ab [□ + m 2 + £R(x)] + 'd-c abcd 4> c (x)(f> d (x ) 


A 


abed ' 


8 4 (x — x') 


V^g 7 


(5.33) 


Similarly, we obtain the mean-held equation for 0 at leading order in the 1/N expansion, 

1 


□ + m 2 + + ^0 2 + ^ G ++ (x , x) ) (f)(x ) + O 


IV 


= 0 , 


(5.34) 


where we note that G ++ (a;, x) = G a b(x,x) for all a,b , which can be seen from Eq. ( |5.32|) ; 
therefore, to get a consistent set of dynamical equations, we need only consider the ++ 
component of Eq. (|5.32| ). It should also be noted that G ab (x,x ) is formally divergent. 
Regularization of the coincidence limit of the two-point function and the energy-momentum 
tensor is necessary. Multiplying Eq. ( |5.32| ) by G and integrating over spacetime, we obtain 
a differential equation for the ++ CTP Green function, 


A - hX 

□x + m 2 + ^R(x) + -4> 2 {x) + — G ++ {x, x) ) G ++ (x , x') + O 


1) = 6\x - 

(5.35) 


8 Note that the index b is not to be summed in the right-hand side of Eq. ( [5.3C ), and the c 
subscript on □ and A is a CTP index. 
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where boundary conditions must be specified on G 


++• 


Equations ( 5.34 ) and ( 5.35 ) are the covariant evolution equations for the mean field 0 and 
the two-point function C++ at leading order in the 1/IV expansion. Following Eq. (|5.13| ), we 
denote the coincidence limit of hG ++ (x,x) by (<^|). With the inclusion of the semiclassical 
gravity field equation i [l. 1 l| j. these equations form a consistent, closed set of dynamical 
equations for the mean field 0, the time-ordered fluctuation-field Green function G ++ , and 
the metric g^ v . 

The one-loop equations for 0 and G can be obtained from the leading-order equations by 
setting h = 0 in Eq. ( |5.35|) , while leaving the mean-field equation (|5.34| ) unchanged. In the 
Hartree approximation, the gap equation is unchanged from Eq. (|5.32|) , and the mean-field 
equation is obtained from Eq. ( 5.34 ) by changing h —> 3 h [^]. The principal limitation of 
the leading-order large -N approximation is that it neglects the setting-sun diagram which is 
the lowest-order contribution to collisional thermalization of the system [|34[]. The system, 
therefore, does not thermalize at leading order in 1/IV, and the approximation breaks down 
on a time scale r 2 which is on the order of the mean free time for binary scattering 075 


(collisional thermalization processes in reheating the post-inflationary Universe is discussed 
in a subsequent paper j50H). 

Let us now use Eq. ( |4.13| a) to derive the bare semiclassical Einstein equation for the O(IV) 
theory at leading order in 1/IV. This equation contains two parts 5S G /5g^ and 5T/5g^. 
The latter part is related to the bare energy-momentum tensor (T0,) by Eq. ( |4.15| ). At 
leading order in 1/IV, (T^ u ) is given by a sum of classical and quantum parts, 


\N 

<r„„} = r« +r« <rf)V> 


(5.36) 


where we define the classical part of (T0,) by 


T c = N 

flV 


(1 - 2- g J 9 P v9 pa ^-A-,a ~ 2£0 ;/il ,0 

+2^000 - £G^0 2 + ^g pu (m 2 + j4> 2 j 0 2 


(5.37) 


and the quantum part of (T0,) by 


T® = Nk lim 

LIU . 

X'^X 


(1 - 20V M V/ + ( 20 - - ) g^g pu ^ P V a - + 2 


G ++ (x, x') > + 0(1). (5.38) 


+ -g jlv ^ m 2 + — 0 2 + — G ++ (x, x 1 ) 

The above expression for T G W is divergent in four spacetime dimensions, and needs to be 
regularized or renormalized. The energy-momentum tensor in the one-loop approximation 
is obtained by neglecting the 0{h 2 ) terms in Eq. (|5.38|) . It can be shown using (|5.35|) that 
the energy-momentum tensor at leading order in the 1/IV expansion is covariantly conserved, 
up to terms of order 0(1) (next-to-leading-order). The bare semiclassical Einstein equation 
is then given (in terms of (T^ u ) shown above) by Eq. ([4.14| ). 

At this point we formally set IV = 1 since we are not including next-to-leading-order 
diagrams in the 1/IV expansion. This can be envisioned as a simple rescaling of the Planck 
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mass by \Z~N, since the matter field effective action T is entirely O(N). We now turn to the 
issue of renormalization. 


D. Renormalization 

To renormalize the leading-order, large-fV, CTP effective action in a general curved 
spacetime, one can use dimensional regularization |67|], which requires formulating effective 
action in n spacetime dimensions. This necessitates the introduction of a length parameter 
/i -1 into the classical action, A —> A/i 4 ^, in order for the classical action to have consistent 
units. As above, we maintain the restriction g+' = g^ u = g^, and we suppress indices inside 
functional arguments. 

Making a substitution of the gap equation into the leading-order, large- N, 2PI effective 
action, we obtain 


r[<MT = s F [0,<H + ^triii (g- 1 ) 


h 2 NXy 4 ~ n 


IM 


d n xy/^gc abcd [ G a b(x , x)G cd (x, x)], (5.39) 


in terms of the CTP propagator G a b(x,x') which satisfies the gap equation 


(G~Y = * {V x c ab + X ab (x)) S(x - x ') 

in terms of a four-component “effective mass” 

A/i 4 n 




X ab (x) = (m 2 + £R) c ab H--— c abcd [4) c 4) d + hGcd(x, x)]. 


(5.40) 


(5.41) 


The divergences in the effective action can be made explicit with the use of the heat kernel 
K a b (x,y,s ) [|68| , |67| , ^9i . Let us define K a b (x,y,s) which satisfies 


(5.42) 


— 7 V ‘ + [ 4" - v —G ')“(!, z)K\(z, y, s) = 0, 


’ M 


with boundary conditions 


K a b (x,y-0) = 8 a b S(x — y) 


V 9y 


(5.43) 


at s = 0 H59fl . From Eqs. ( |5-42|) and (|5.39|) it follows that K + _ = K~ + = 0 for all x, y , and 
s, and that K + + (K~_) is a functional of 0 + (0_) only. The CTP effective action can then 
be expressed as 


P[0, <r] = r to[k, 9 n - r^-, (5.44) 

in terms of a functional r io on M defined by 
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r+„[0 + ,<n = s'P+.sT] 


ihN 


ds 


d n xyf~g / — K + + (x, x] s) 


* JM 

h 2 N\n A ~ n 


d"Xy/=7) 


IM 


UO 


dsK + + (x, x ; s) 


(5.45) 


and similarly for r io . It follows from Eq. (|5.42|) that K + + (x,x] s)[(p + ] is exactly the same 
functional of <p + as K~_(x, x; s) [</>_] is of </>_; we denote it by K(x,x;s)[(f)], where 0 is a 
function on M. 

The divergences in the effective action arise in the small-s part of the integrations, so 
that in the equation 



ds 

s 


K(x, x ; s ) 



K(x, x ; s ) 


+ 



ds 

s 


K(x, x ; s) 


(5.46) 


only the first term on the right-hand side is divergent. Using the s —> 0 + asymptotic 
expansion for K(x,x;s ) [^9[, one has (for a scalar field, such as the unbroken symmetry, 
large- N limit of the O (N) model) 


K(x, x ; s ) rs./ (4tts) 2 ^ 


s m a n 


w, 


(5.47) 


m =0 


where the a n {x) are the well-known “Hamidew coefficients” made up of scalar invariants of 
the spacetime curvature [jl],|68|]. The divergences then show up as poles in l/(n — 4) after 
the s integrations are performed. They have been evaluated for the A<h 4 theory in a general 
spacetime by many authors (see, e.g., lp0| . [70| . [7T| ) and in the large- N limit of the O (N) model 


0. At leading order in the 1/N expansion, the renormalization of A, £, m, G , A, b , and c 
is required, but no held amplitude renormalization is required. 
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FIGURES 



FIG. 1. Diagrammatic expansion for Ti. Lines represent the propagator x'), and vertices 

terminating three lines are proportional to cj). Each vertex carries spacetime (x) and CTP (+, —) 
labels. 



FIG. 2. Diagrammatic expansion for T 2 - Lines represent the propagator G, and vertices are 
given by 51 1 . The vertices terminating three lines are proportional to cj). 
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